Network Model for a 2D Disordered Electron System with Spin-Orbit Scattering 
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We introduce a network model to describe two-dimensional disordered electron systems with spin- 
orbit scattering. The network model is denned by a discrete unitary time evolution operator. We es- 
tablish by numerical transfer matrix calculations that the model exhibits a localization-delocalization 
transition. We determine the corresponding phase diagram in the parameter space of disorder scat- 
tering strength and spin-orbit scattering strength. Near the critical point we determine by statistical 
analysis a one-parameter scaling function and the critical exponent of the localization length to be 
v = 2.51 ±0.18. Based on a conformal mapping we also calculate the scaling exponent of the typical 
local density of states olq — 2.174 ± 0.003. 

PACS numbers: 71.30.+h; 73.23.-b; 71.70.Ej; 64.60.-i 



i 

-a 
c 

o 
o 



> 

(N 

m 
m 
O 
oo 
On 

-i— > 
c3 



i 

C 

O 

o 



X 



I. INTRODUCTION 

Recently, localization-delocalization (LD) transitions 
in 2D disordered electron systems in the absence of a 
magnetic field were observed by several groups [0-0- 
These results are in contrast with the scaling theory 
for non- interacting electrons |^|, which predicts that all 
states are localized in two dimensions and in the absence 
of spin-orbit interaction (SOI). Now, a new discussion 
has started on this topic with the emphasis on the effects 
of electron-electron interaction and spin-orbit interaction 

HI- 

It is known that both types of interactions could be re- 
sponsible for the existence of a LD transition. In the case 
of SOI, general arguments H] and perturbation theoret- 
ical calculations in the weakly disordered regime [ p|fl6| 
yield a positive correction to the conductance. This 
quantum interference effect requiring time reversal invari- 
ance is known as weak anti-localization. In the present 
work we focus on the detailed examination of a 2D non- 
interacting electron system with SOI. For these purposes 
we formulate a scattering theoretical network model for 
such a system. 

In a recent paper [|l7j is was shown that scattering the- 
oretical network models (NWMs) are well suited to de- 
scribe mesoscopic disordered electron system. In general 
such a NWM can represent any system of coherent waves 
propagating through disordered media. It consists of a 
network of unitary scatterers connected by bonds. The 
arrangement of scatterers and bonds defines the topology 
of the NWM, which can be described by a connectivity 
matrix. In our work we have chosen a simple case, where 
the scatterers are located on the sites of a quadratic grid, 
so each of them has four nearest neighbors. Each bond 
consists of 2n links, n for each direction, where n = 1 for 
waves without and n > 1 for waves with internal degrees 
of freedom (cf. Fig. |l|). In the case of electron waves 
a complex number is attached to each link representing 
the probability amplitude at this position. The set of all 
amplitudes defines the quantum mechanical state ty(t) 



at time t. One step of time evolution is then given by a 
unitary operator hi, 



#(t + l) =W*(i). 



(1) 



This time evolution operator is determined by all the 
scatterers in the NWM. Each scatterer maps An incoming 
channels to An outgoing channels conserving the current 
and is therefore represented by a unitary An x 4n-matrix. 
The disorder is in general simulated in two ways: first 
by multiplying the amplitude on each link with a com- 
plex random phase factor e 1 ^ with <f> randomly chosen 
from [0, 2-7r[ simulating the random distances between the 
scatterers and secondly by taking random values for the 
parameters that parameterize the matrix representation 
of the scatterers simulating the random strengths of the 
scatterers. Of course, both random choices have to be 
compatible with the symmetry properties of the system. 
We distinguish 2D electron systems with time reversal 
symmetry (02NC) and without time reversal symmetry 
(U2NC), both without spin degrees of freedom, and sys- 
tems with time reversal symmetry and spin degrees of 
freedom (S2NC). The "2" refers to the space dimension, 
the "O", "U" and "S" mean " orthogonal" , " unitary" and 
" symplectic" , which refers to the corresponding univer- 
sality classes of random matrix theory and the letters 
"NC" indicate that all these systems are " non-chiral" , 
which means that no orientation is preferred as would be 
in presence of a strong magnetic field. 

The two former models have been examined exten- 
sively in fL8| . In this work a reflection, a transmission 
and a deflection coefficient were introduced, which pa- 
rameterize the scattering matrices. Furthermore an elas- 
tic mean free path was defined in terms of these coeffi- 
cients. It was concluded that all states are localized in 
the 02NC-/U2NC-NWM. 

In the present work we investigate the S2NC-NWM. 
We find a parameterization for the matrix representa- 
tion of the spin scatterers and introduce a spin scatter- 
ing strength. The calculation of the localization length 
by the transfer matrix method allows us to detect the 
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FIG. 1. Topology of a general network model. Squares are 
scatterers and lines are bonds. 



LD transition and to determine the scaling function. In 
order to quantify the scaling exponent v of the correla- 
tion length we use a fit procedure Jl9| . We fit the scaling 
function and the critical exponent in two steps respecting 
the correlations of the data. Additionally, we apply a \ 2 - 
test to estimate the confidence of the fits. Determining 
the critical value of the localization length we find the 
scaling exponent ao of the typical local density of states 
using a conformal mapping |2Cf |. 

This paper is organized as follows: In Sec. |l] we in- 
troduce the network model by explicitly constructing the 
scattering matrices. Sec. [II contains the transformation 



to the transfer matrices and summarizes general aspects 
of LD transitions. A detailed description of the methods 
of data evaluation forms the content of Sec. fV|. The 
discussion of the results is presented in Sec. |V followed 
by a short summary in Sec. [VT| . 



FIG. 2. Topology of the network. The potential scatterers 
(white) change the direction, the spin scatterers (grey) the 
spin of the electrons. 
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FIG. 3. Potential scatterer: The eight incoming channels 
7f are mapped to the eight outgoing channels Of . 



II. NETWORK MODEL 
A. Topology 

There are two different types of scatterers in the S2NC- 
NWM: potential scatterers (PSs) changing only the elec- 
tron's direction and spin scatterers (SSs) changing only 
the electron's spin. The network consists of a regular 2D 
quadratic grid of potential scatterers each of them con- 
nected to the four next neighbors by bonds. On each 
bond a spin scatterer is placed leaving the electron's di- 
rection unchanged (cf. Fig. 

B. Potential scatterers 

There are four channels or links within each bond, two 
of them for incoming and two of them for outgoing states 
with spin up and spin down, respectively. The electronic 
state is represented by complex numbers (amplitudes) on 
each link. Consequently, each PS maps eight incoming 
channels I? to eight outgoing channels Of (a G {+,—}, 



i G {1,2,3,4}) and thus can be represented by a 8 x 8- 
matrix S po t- With the definition of the geometrical ar- 
rangement of the channels shown in Fig. ^ this mapping 
is defined as follows: 
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Due to current conservation, 

Ei^i^Ei^i 2 ' ( 3 ) 

each scattering matrix has to be unitary, 

S P ot • S^t = Is, (4) 

where Is denotes the 8x8 identity matrix. Additionally, 

the scatterers are time reversal invariant. Both proper- 
ties yield the matrix to be symmetric, 

Spot = Spot ; (5) 
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FIG. 4. Definition of the reflection coefficient r, the trans- 
mission coefficient t and the deflection coefficient d. 



where T denotes the transpose. 

For convenience we choose the potential scatterers to 
be isotropic, i.e. they are invariant under rotations by 
multiple angles of tt/2. With these restrictions each scat- 
tering matrix S po t can be parameterized in the following 
way pq|: 



->pot 



with 



Jpot 



and 
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(6) 



(7) 
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(8) 



Here 11 2 denotes the 2x2 identity matrix and ® is the ten- 
sor product. The real parameters r, t, d denote the reflec- 
tion, transmission and deflection (right and left scatter- 
ing) coefficient, respectively (cf. Fig. [|). If we choose r 
and t as independent parameters for S po t , the real phases 
<f> r , (fit and the deflection coefficient d are related to them 
due to unitarity and time reversal symmetry, 



2\d\ 2 + \t\ 2 =1, 



\r\ ■ \d\ cos (fi r 
\r\ ■ \t\ cos((fi r - (fit) 



= -\t\ 
2 



|d| cos (fit 



(9a) 
(9b) 
(9c) 



Furthermore, two restrictions follow from these equa- 
tions, 



r 2 + t 2 



r +t 



< 1, 
> 1. 



(10a) 
(10b) 



The four real phases <fi\ , . . . , <f>4 which are randomly 
chosen from the interval [0,27r[ model the spatial disor- 
der. They can be interpreted as phase factors e 1 ^ for 
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FIG. 5. Spin scatterer: The four incoming channels I" , I" 
are mapped to the eight outgoing channels O a , O" . 



freely propagating electron waves. Consequently, there 
are six independent parameters. But only r and t gov- 
ern the macroscopic properties of the system. For con- 
venience we choose them to be equal for all PSs in the 
network, whereas the phases are randomly taken for each 
scatterer. 



C. Spin scatterers 

The spin scatterers located between the potential scat- 
terers have two incoming and two outgoing channels on 
the left and on the right, respectively, as is shown in Fig. 
||. They can be represented by 4 x 4-matrices S sp : 



O 



S J 



with I 



1+ \ 










' o~ 


1+ 








\ 0- 



Unitarity and time reversal symmetry result in 

1T C T 



and 



S sp = D T ^S- 1 ifD = D T (S s - 1 ) it D, 



(11) 



(12) 



(13) 



respectively. Here the asterisk denotes complex conjuga- 
tion and KD is the time- reversal operator with complex 
conjugation operator K and 
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s 




= 








(14) 



where T2 is one of the basis quaternions [tq,t) 
(to,t 1 ,t 2 ,t 3 ) given by 



tq = I2 and 



t = —\<j 



with Pauli matrices cr — (a x , o~ y , o~ z ). The matrix 



J = 



1 

-1 



(15) 



(16) 



is the symplectic unit matrix. The symmetries ( p_2|) and 
( |T3| ) suggest the following parameterization of the spin 
scattering matrix, 



e llp q 
e iip q 



(17) 
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where the real coefficients qi are restricted by X)fe=o 9fe = 
1 due to unitarity and q denotes the quaternion conju- 
gation ]2ff| . Thus, three independent parameters remain 
which are randomly and homogeneously taken from the 
unit sphere for each SS. The phase ip is randomly taken 
from [0, 2tt\. In random matrix theory the symmetries 
( |l2| ) and (|l3|) correspond to an ensemble of quaternion 
real Hamiltonians, which can be diagonalized by sym- 
plectic transformations. Because of this, we also refer to 
this symmetry as symplcctic symmetry. 

We now define the spin scattering strength by 



= V 1 - qo = 



q'l + qi + ql 



(19) 



This quantity takes values in the interval [0, 1], where s = 
means no spin scattering and s = 1 full spin scattering, 
resulting in full spin relaxation after one scattering event. 
The parameter s is fixed for the whole network. 



D. Parameter space 

In conclusion there are three independent quantities 
building up the three dimensional parameter space (or 
phase space) of the S2NC-NWM. With the restrictions 
( |i"o| ) and s G [0, 1] the possible values (r, t, s) are located 
in a certain volume in R 3 . Fig. |^ shows a cross-section of 
this volume at some fixed value of s. The grey area con- 
tains the allowed values. If the time reversal symmetry 
is omitted, i.e. switching to the U2NC-NWM, the phase 
space (at some fixed s) is the entire quarter of the circle. 

There are three exceptional points in the phase space: 
For r — and t = 1 we have the delocalized fixed point, 
where the electron waves propagate freely. If we take 



r = 1 and t = 0, we have the localized fixed point, where 
transport stops. For r = t = the network is at the 
Chalker-Coddington fixed point |22|, where only left and 
right scattering exists. On the line r 2 + t 2 — 1 the de- 
flection coefficient is zero. Thus, the system splits in 
independent ID subsystems which always show localiza- 
tion. All of these properties are independent of the value 
of s. 

In an elastic mean free path is defined by 



2 r 2 + d 2 



(20) 



The factor 1/2 is a consequence of the diagonal arrange- 
ment, which is shown in Fig. [?]. The unit of l e is a lattice 
constant. Analogously, we define a spin scattering length 
by 



^so := 



1 1 
2~ 



(21) 



This length scale takes values from for maximal spin 
scattering to oo in the case of the absence of spin scat- 
tering. 



III. FINITE-SIZE SCALING 

A. Transfer Matrix Method 

In order to investigate the scaling behavior of the mod- 
eled system we need some scaling variable. Although the 
conductance is the natural choice for a scaling variable in 
this context, the renormalized localization length (RLL) 
A = £/M is, as a self-averaging quantity, much more con- 
venient. Here £ is the localization length of a quasi- ID 



4 



system of width M and can be calculated by the transfer 
matrix method 23,24}], This method yields a sequence of 
Lyapunov exponents in decreasing order, where (due to 
conservation of current density) each value has a partner 
with opposite sign. Beyond this, in presence of SOI each 
value appears twice, because of time reversal symmetry 
(Kramers degeneracy). The smallest positive Lyapunov 
exponent determines the quasi-lD localization length. It 
is always finite due to the finite width of the system. 

In order to be able to apply the transfer matrix method 
we have to convert our scattering matrices S pot and S sp 
into the corresponding transfer matrices T pot and T sp , 
which map channels on the left to channels on the right. 
In contrast to the scattering matrices the transfer ma- 
trices are multiplicative, which means the following: The 
total system is divided into a sequence of elementary sub- 
systems each of them corresponding to a single transfer 
step. The transfer matrix of the total system is then 
given by the product of the transfer matrices of the sub- 
systems ("strip transfer matrices"). 



B. Transfer Matrices in the Network Model 
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(25) 



(26) 



(27) 



respectively. On the other hand, time reversal invariance 
requires 
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pot 



D 
D 



T* 



D 

pot 1 q 



and 



T — D T T* D 



(28) 



(29) 



with D as given in ( |l4| ) . A parameterization of T pot com- 
patible with these restrictions is given by [jig] 



With the channel orientation of Fig. || and Fig. || the 
transfer matrices of the S2NC-NWM are defined by 
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(22) 
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(23) 



respectively. From this definition a diagonal arrange- 
ment of the whole network results, as is shown in Fig. 
|?]. The arrows at the top and the bottom of the figure 
indicate the periodic boundary conditions which we have 
chosen. The natural width unit in this arrangement is 
a pair of diagonally neighbored transfer matrices, which 
corresponds to a channel number of A c = 8. Therefore 
the bold printed part of the picture represents a strip 
transfer matrix of width M — 2, channel number iV c = 16 
and unit length L = 1 (horizontal direction). 

In the language of transfer matrices the conservation 
of current density writes as pseudo-unitarity: 



T V TT 

1 pot^z 1 po t 



(24) 



with 



T 



pot 



/ a* 7 (3* S 
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(30) 



with 



d 

A? 

(re 1 *- - te ict>t )d 



te 1 ' 



A 



5 = 



d 2 - r&^te 1 ' 
A 



(31) 



A = d 2 - (te^*) 2 



For T sp we find 



T s 



e iv q 
e-'fq 



(32) 



with q as in Eq. (18). We omitted the four phase factors 
in T po t because they can be combined with those of T sp 
to a resulting phase factor ip in T sp . 



C. Localization-Delocalization Transition 

The scaling behavior of the RLL A determines whether 
the system is localized or delocalized. If A shrinks with 
increasing system width M the system behaves like an 
insulator, if it grows, the system is metallic. At the LD 
transition the RLL is independent of the system width. 
To ensure that A is in fact a scaling variable one has to 
find a scaling function that is a function of only the ratio 
of the correlation length £ c and M, 



5 



A(M) = / 



M 



or logarithmically 



lnA(M) = / lnM-m£, 



(33) 



(34) 



Equivalent to the existence of a scaling function is the 
formulation of a flow equation with a /3-function that is 
a function of In A only, 



/3(lnA) 



din A 
d In M ' 



(35) 



Fig. H shows a qualitative picture of the /3-function in dif- 
ferent dimensions. Due to Ohm's law the limiting value 
of /3 for large conductance (or A) is d — 2. Without SOI 
there is a negative correction to the conductance caused 
by weak localization |^f|. The /3-function is always neg- 
ative in 2D and therefore all states are localized. In 
the presence of SOI the correction to the metallic con- 
ductance changes the sign due to weak anti-localization 
E p"6 |. Consequently, in 2D and in the large conduc- 
tance imit the /3-function converges to zero from above. 
Since in the strongly disordered regime all states are ex- 
ponentially localized, the existence of a LD transition in 
the 2D symplectic case follows from simple scaling argu- 
ments. At the critical point, where the /3-function is zero, 
the correlation length shows power low scaling with the 
critical exponent v, 



Ur)=£\r-r*\- 



(36) 



where r is a system parameter, e.g. the reflection co- 
efficient of the NWM, and r* is its critical value. Here 
the correlation length is the fictitious system width up to 
which the system is in the critical regime. In the local- 
ized regime this is just the quasi-lD localization length 
obtained by the transfer matrix method in the thermo- 
dynamic limit, 



£ c = lim £(M)=£ 

M— »oo 



(37) 



In EQ] it was shown that the typical local density of 
states (LDOS) is an appropriate choice for the order pa- 
rameter of the LD transition. 



Ptyp 



e (ln P (B,r)> K 



(38) 



with p(r) — \ip(E, r)\ /A(E), energy E and local level 
spacing A(_E). (...) denotes disorder average and (3 p is 
the critical exponent of the order parameter. Further- 
more, ptyp shows power law scaling at the critical point 
with exponent d — olq, 



Ptyp oc L 



d-a 



(39) 



where L d is the volume of a d-dimensional cube and ao 
is a scaling exponent, which is known from multifrac- 
tal analysis of critical wave functions. In 2D, i.e. for a 




d=2 in A 



FIG. 8. Qualitative picture of the /3-function for d = 1, 2, 3. 
In the 2D case there can be a transition due to weak 
anti- localization. 



square system, this scaling exponent is linked to the crit- 
ical value of the quasi- ID RLL by a conformal mapping 
argument BQ], 



A* 



7r(a - 2) 



(40) 



With the knowledge of v and ao the critical exponent of 
the LDOS is given by 

(3 p = v(a - 2). (41) 



IV. METHODS OF EVALUATION 

A. Fit Procedure for the Scaling Function 

In our work we follow the method introduced in [ jl9| , 
which not only fits the scaling function to the data of the 
RLLs but also uses a x 2 test to check the confidence of 
the fit. 



According to (34) we want to fit the scaling function 
/ to the logarithms of the RLLs, which depend on the 
njvf system widths {Mi, . . . , M nM } and the n r system 
parameters {ri, . . . , r Ur }. So we have n\ — hm • n r data 
points {Ai , . . . , A„ A }. For abbreviation we now introduce 
the following vectors, 



Y := 



and 



/ In Ai 



V lnA„ A 



lnMi-]n£ c (n) 



X := 



lnM nM - ln£ c (r„J 



AlnAi 



AlnA r 



(42) 



(43) 



the latter being the vector containing the errors of the av- 
erage obtained by the transfer matrix method. Assuming 



G 



the data to be statistically independent the correspond- 
ing correlation matrix Ca given by 



(C A )„ := (Ei-Ej) 



(44) 



is diagonal. Here (. . .} denotes the mean value. 

Following the procedure in Ref. jl9| we make an ansatz 
for the scaling function by a linear combination of Cheby- 
shev polynomials, 



F(xi) 



2 



N 

E 

k=l 



CkTk(Xi 



(45) 



which gives a polynomial of degree N. Omitting the 
tilde over the Xi indicates the argument being rescaled 
to [—1, 1]. In this interval the Chebyshev polynomials are 
orthogonal and have simple behavior at the edges. 

The smallest of the parameters ln£ c is fixed by re- 
quiring ln£ c (n) = for the delocalized branch and 
ln£ c (r„ r ) = for the localized branch, respectively. So 
we have the tiq = n r + N parameters 

/ InCc(ri) \ 







for the localized and 



:= 



hi^ c (r„ r _i) 
Co 

V c « J 

( InCcM \ 

hi£ c (r„ r ) 

CO 



(46) 



V 



c/v 



(47) 



/ 



for the delocalized branch, respectively. These parame- 
ters have to be fitted with respect to the data Y , Hence, 
the n\ values of the fit function can be written as 



F(X;®) 



F(X i; &) 



F(X nA ;@) 



(48) 



We use the method of the least-squares fit, i.e. we 
have to minimize the sum Se of the weighted quadratic 
deviations, which means solving 



dS e 
d® 







with 



S @ = (r-F(X;0)) T CX 1 (^-F(X;0) 



(49) 



(50) 



Since F is non-linear in the parameters this proce- 
dure leads to a system of ne coupled non-linear equa- 
tions. Therefore So is minimized directly by a numerical 
method. Starting with some estimated initial values for 
the logarithms of the correlation lengths we successively 
optimize the c,; and the ln£ c (7-j). If the data are compati- 
ble we will have convergence, and thus we can stop when 
a chosen accuracy is reached. The result then is a set of 
coefficients c, which defines the fit function and a set of 
optimized values for the correlation lengths. 

According to Ref. |2(| the correlation matrix of the 
parameters is given by 



C<-> 



(51) 



where Fe is the Jacobian of F with respect to 0, 

dF 

Fe := (52) 



Usually, as errors of the parameters one takes the diago- 
nal elements of the error matrix E, which is defined by 



S'e 



n\ - ne 



-C 



(53) 



B. Testing the Fit of the Scaling Function 

A converging fit procedure doesn't guarantee that the 
errors of the numerical data are actually compatible with 
the obtained fit function. Therefore, in addition to this fit 
procedure we apply a % 2 test to estimate the confidence of 
the fit. We make the essential assumption that the data 
Hi are normally distributed about f(xi] ©) with variances 
ej 2 . Consequently the quantity S'e has to be distributed 
as x 2 with ha — n® degrees of freedom. This distribution 
has the estimated value n& — ne and the variance 2(?ia — 
ne). A suitable measure for the confidence of the fit then 
is the normalized deviation of S from the estimated value 



■S'e - (»a - ne) 
v /2(n A - n e ) 



(54) 



If |Ae| < 1 it is safe to assume that the fit is trustworthy. 
But if | A© | takes values much larger than 1 it is very 
unlikely that the data Yi are normally distributed about 
F(Xi;Q) which indicates systematic errors. In this case 
the fit has failed, we have to give up the assumption of 
one-parameter scaling and further calculations, e.g. of 
the critical exponent, don't make much sense. 

It should be noticed that rescaling of the variance ma- 
trix S by a factor b, S = b ■ S, results in the reciprocal 
rescaling of Se, Sq — b~ 1 Se. Since Ae depends sensi- 
tively on S'e, especially if n\ is small, one should carefully 
consider, how to determine the errors of the raw data. 
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C. Fitting the Critical Exponent v and r* 
and Testing the Fit 

In order to determine the critical exponent of the cor- 
relation length we apply the same idea as before with 
respect to Eq. ((3^) . Particularly, we simultaneously deal 
with both branches of the scaling function, i.e. we as- 
sume the critical value r* and the critical exponent v to 
be the same in the localized and the delocalized regime. 
Only the prefactor can take different values, which will 
be denoted as £c' loc an d £°> deloc for the localized and the 
delocalized regime, respectively. Taking the logarithm 
Eq. (pq) writes as 



ln&(r<) = lng> l ° c -vln\n-r 



and 



InCcO*,) = ln£ 



0,dcloc 



ybxl 



(55a) 



(55b) 



for Ti in the localized and delocalized regime, respectively. 

The n r values ln£ c (ri) are the results of the foregoing 
optimization, the arguments are the values In |r*j — r* \ 
and the four parameters that have to be optimized are 
£c' 1oC j £°' dcloc , v and r*. Introducing the vectors 



(56a) 



/ m£c(ri) \ 




f In |n 




, x := 




V lnCc(r„J / 




\ m \ r n, 



( ln£c°' loc \ 

m ^0,dcloc 
V 

r* 



m £Uoc 

and 9 := | ln£?< doloc 
v 



(56b) 



and comparing them with Eq. (J55]) the fit function can 
be written as 



with 



W 



f(x;9)=\N(x)0, 

( 1 -hx\n-r*\ \ 

10 — In \r nr loc — r* | 
1 -In |r„ r loc+ i - r*\ 



(57) 



\o 1 



- In |r„ 



(58) 



1 J 



Here n rj i oc is the number of values which belong to the 
localized regime. 

The correlation matrix C^ c of the data yi is the upper 
left n r x n r submatrix of Ce (Eq. (pi])) obtained by the 
fit of the scaling function. Thus, the sum of the quadratic 
deviations is 



S e = 



y-f(x-9)Yq c 1 {y-f(x-9)y 



(59) 



Since the fit function is a linear function in the argument 
x, we can analytically solve the minimization problem 



dSg 

89 



0. 



which leads to 



9= (\N T C7 1 \N) W T C7V 



(60) 



(61) 



In contrast to that the parameter r* has to be optimized 
numerically since it appears non-linear ly in Eq. j5^). 
Giving some starting value for r* one iteration step con- 
sists of successively optimizing 9 and r* . 

The 4x4 correlation matrix of the four parameters is 
given by 



Co = (FjC^F 



where 



F • 9f 

Finally we get the error matrix 

F Se r 
te = -L 

n r — 4 



As in Eq. (54) we use the quantity 
S e - K - 4) 



V 2 K - 4) 



(62) 



(63) 



(64) 



(65) 



to test the confidence of the fit. |A#| should take val- 
ues of about 1 or smaller to verify the assumption that 
the values yi are normally distributed about the fit func- 
tion, which is a straight line in this case. Note that it 
is important to take the correlations between the ^ c (^i) 
introduced by the previous fit procedure into account in 
the present analysis. Thus, a simple linear regression will 
not give the correct results. 



D. Determination of A* 



The fit procedures introduced in Sects. IV A and IV C 
are not the most direct way to determine v from the 
raw data. Instead, one can fit the RLLs as functions of 
r — r*| with a width dependent scaling factor M x l v 



A(M; r) = h (m 1 ^ \ 



(66) 



which is a consequence of Eqs. (|33|) and (|36|). This pro- 
cedure leads to a continuous curve, because there is no 
splitting in two branches as in the logarithmic case. Fit- 
ting the scaling function in this manner allows a direct 
evaluation of A*. Since at r = r* the argument of h is 
zero, one only has to calculate the function value: 
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S 0.30 - 

< 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 9. Renormalized localization length for strong disor- 
der, r — 0.8, t = 0.4 which corresponds to l c ~ 0.35, depend- 
ing on the spin scattering strength s and the the system width 
M. 



FIG. 10. Renormalized localization length for r = 0.55, 
t = 0.6 which corresponds to l c ~ 1.1, depending on the spin 
scattering strength s and the the system width M. 



A* = h(0). (67) 

Following the law of propagation of errors we get for the 
error 

AA* = AA| r=r » = M^Ah'iOy*, (68) 

where hi denotes the first derivative of h with respect to 
the argument. 

V. RESULTS AND DISCUSSIONS 
A. Localization Lengths 

We calculated quasi- ID localization lengths for system 
lengths up to L = 2 ■ 10 5 and widths from M — 2 up 
to M = 32, which corresponds to channel numbers from 
N c = 16 to N c = 256. The corresponding errors AA 
are the errors of the average, which vanish for L — > oo. 
We have chosen such system lengths L that the relative 
errors AA/A take values of about 0.1 % to 1 %. 

Fig. [| shows the renormalized localization lengths for 
r = 0.8 and t — 0.4 decreasing with increasing system 
width for the whole range of the spin scattering strength. 
So the system is in the deeply localized regime. This 
matches with the fact that the mean free path (see Eq. 
(pp|)) takes a value of l c ~ 0.35 in units of the lattice 
constant. So the reflection is too strong to allow extended 
states at all. 

Contrary to this, the next example exhibits a LD tran- 
sition at s ~ 0.3, which can be seen in Fig. |l^. With 
r = 0.55 and t = 0.6 corresponding to l c ~ 1.1 the ef- 
fect of weak anti-localization causes the existence of ex- 
tended states, if s is strong enough. The intersection of 
the curves clearly indicates the LD transition. However, 



their slope is rather small compared to the error bars pre- 
venting an accurate scaling analysis close to the critical 
point. 

The third example shown in Fig. [ll] belongs to the 
values t = 0.8 and s = 0.4, where r varies from r — 0.48 
to r = 0.56. It exhibits another kind of problem. Al- 
though curves for different system widths M intersect, 
the points of intersection systematically depend on the 
width. The larger M is the closer are the points of in- 
tersection for curves of neighboring values of M. It is 
obvious that there exists a limiting point, which would 
be the true critical point. The observed deviations are 
the consequence of finite-size effects, which vanish in the 
thermodynamic limit. By investigating A at a lot of dif- 
ferent points in parameter space we have seen that the 
deviations due to finite-size effects are larger for more de- 
localized systems. Actually, there is only a small area in 
parameter space that is suitable for a quantitative anal- 
ysis of the LD transition. 

B. Phase Diagram 

We determined a phase diagram for the LD transition 
by using the scaling behavior of the RLL. More precisely, 
we calculated A(Mx = 4) and A(M 2 = 8) for a lot of 
pairs (r,i) with s = 0.01, 0.02, 0.05, 0.1, 0.4 and 1. In 
order to get the critical line in the (r, t)-subspace with 
fixed s we decided the point (r, t, s) to be in the localized 
and delocalized regime, if 

A(Mx) - AA(Mi) > A(M 2 ) + AA(M 2 ) (69) 

and 

A(M 2 )-AA(M 2 ) > A(Afi) + AA(Mi), (70) 
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FIG. 11. Renormalized localization length in the vicinity 
of the localization-delocalization transition, t = 0.8, s = 0.6 
which corresponds to l e ~ 2, depending on the reflection co- 
efficient r and the the system width M. The arrows mark the 
section of curves corresponding to neighbored values of M. 



respectively. In the case that both conditions failed, we 
considered the point in the parameter space to be criti- 
cal. By that procedure we got a critical region, i.e. the 
separating line had some finite width. Fig. ^ shows 
the resulting phase diagrams for some intersections of 
the parameter space at the above declared fixed values 
of s. The white area marks the localized the grey one the 
delocalized phase. The drawn critical line is the optical 
interpolated center line of the critical region. 

The region of the metallic phase shrinks with decreas- 
ing spin scattering strength. This is due to the fact that 
the weak anti-localization then becomes less effective in 
preventing Anderson localization. At s = the system 
changes universality class from symplectic to orthogonal 
symmetry, a fact that could be verified by comparing the 
values of the localization lengths with those in Ref. |L8] . 
On the other hand even a very small value of s gives rise 
to a certain delocalized phase, if r is small and t is large 
enough. Of course, the larger l e , i.e. the larger t and 
smaller r, the more easily extended states occur. But 
even in the presence of full spin scattering only about 
the half of the area of the parameter space belongs to 
the metallic phase. This is due to the fact that parameter 
values (r, t) belonging to the localized phase correspond 
to too strong disorder resulting in too strong localization 
to be broken by weak anti-localization. 

It should be noticed, that for t > 0.6 the shape of the 
phase boundary is influenced by finite-size effects. So the 
phase diagram can only serve as a qualitative picture of 
the LD transition. In order to improve on the phase dia- 
gram, one has to consider much larger systems, which is 
very computer time consuming taking such a large num- 
ber of data points into account. Nevertheless, the lower 
part of the boundary, i.e. the region close to the line 
r = t (dotted line in the figure), is suitable for quantita- 





FIG. 12. Phase diagrams for the localization-delocalization 
transition for cross-sections at s = 1, 0.4, 0.1, 0.05, 0.02, 0.01. 
The grey area shows the delocalized, the white the localized 
regime. The dotted line corresponds to r — t. 



tive investigations, as will be shown in the following. 



C. Scaling Function 

We determined the scal ing fu nction by the fitting pro- 
cedure described in Sec. IV A for t = 0.6, s = 0.4 and 
r 6 [0.52,0.62]. In this small region of the phase space 
the corresponding curves A(M; r) (see Fig. |l3j) are very 



well suitable for a quantitative analysis, because of their 
strong r dependence and the absence of noticeable finite- 
size effects (l c ~ 1). 

Fig. |l4| shows the scaling function with the upper 
branch belonging to the metallic and the lower branch be- 
longing to the localized regime. The curves represent the 
fitted Chebyshev polynomials. The data points are the 
raw data shifted by the fitted values of £ c (r) . We omitted 
the data with M = 2 in the localized and M = 2 and 
M = 4 in the delocalized regime, because these values 
showed systematic deviations due to finite-size effects. 
Also data that are too close to the critical point were 
omitted. For the remaining data the confidence test of 
the fit gives Ae = —0.30 and Ae = 0.27 for the localized 
and delocalized branch, respectively. So the assumption 
of one-parameter scaling is very well confirmed. 

Fixing the fit parameters ln£ c by setting ln£ c (r = 
0.52) = for the delocalized and ln£ c (r = 0.62) = 
for the localized branch (circles in Fig. |l4|) the proce- 
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0.52 0.54 0.56 '* 0.58 0.60 0.560 0.565 0.570 0.575 0.580 



FIG. 13. Renormalized localization length for t = 0.6, 
s — 0.4 which corresponds to l e ~ 1, depending on the re- 
flection coefficient r and the the system width M. 




— fit function kldoc. I 
Or- 0.52 

□ r=0.53 
+ r = 0.535 
Ar=0.54 
ir- 0.545 
Or- 0.55 
*r- 0.555 

— fit function (loc.) 
* r - 0.59 

Or. 0.595 
xr-0.60 
Ar= 0.605 
+ r-0.61 

□ r=0.615 
Or =0.62 



0.0 1.0 2.0 3.0 

111 M - In £;.(/ ) 



FIG. 14. Scaling function for 

the localization-delocalization transition. The upper branch 
corresponds to the delocalized and the lower to the localized 
phase. 



dure has converged with an accuracy of 0.1% after about 
50 iterations. The starting values have a radius of con- 
vergence of about 5. So a rough estimate is sufficient for 
convergence. 

The Chebyshev polynomials used are of fourth order. 
With a lower order it is impossible to fit the curves (as 
indicated by the figure of merit Ae), whereas with an 
order higher than 6 the fitted curves start to follow the 
fluctuations of the data points, which results in a non- 
physical behavior of the scaling function. With an order 
between four and six there is no significant difference in 
the results. 



D. Critical Exponent v and Critical RLL A* 

In order to determine the critical exponent v of the 
correlation length we used the fit procedure presented 



FIG. 15. Dependence of v on the assumed value of r* 
best estimate is found for r* = 0.571 



The 



in Sec. [IV C| . With a starting value r* 6 [0.53,0.59] 
the procedure converges. After about 10 iterations the 
corrections are smaller than 1%. The results of the fit 
are 

v = 2.51 ±0.18 and r* = 0.571 ± 0.002 . (71) 
The prefactors take the values 

In 4°' loc = -7.55 ±0.42 (72a) 

and 



Q.deloc 



= -7.50- ±0.43. 



(72b) 



The confidence test of the fit gives Ag = —0.47, showing 
its high quality. It is very important to stress, that the 
given errors (Eq. j64[)) are not independent. They have 
to be interpreted considering the correlation matrix 



C„ = 



/ 0.0310 -0.0737 -0.0739 -0.0004 

-0.0737 0.1782 0.1766 0.0009 

-0.0739 0.1766 0.1853 0.0009 

V -0.0004 0.0009 0.0009 5.5 ■ 10~ 6 



(73) 



which shows that the different values are highly corre- 
lated. 

Several reference values for v have been published in 
the last years p7|- 0O| . The most recent calculations yield 
v = 2.75 ± 0.1 lljjr^ = 2.5 ±0.3 §| and v = 2.32 ±0.14 
p3| . Within the errors our value agrees with these. We 
note that v is very sensitive to slight variations of r* . 
This is seen by fitting v for fixed values of r* . As is shown 
in Fig. [l^ v changes by 30% if r* changes by about 3%. 
The difficulties in obtaining a credible value for v were 
the r easons to employ the iterative fit procedure of Sec. 

ivcj. 

In order to determine the critical RLL we used Eqs. 
© and @. The result is 



A* = 1.83 ±0.03. 



(74) 
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Finally, Eq. ( ff"o| ) yields the value for the scaling exponent 
of the typical LDOS 



a = 2.174 ±0.003. 



(75) 



This value agrees well with the result of Schweitzer [|34| 
ao = 2.19 ± 0.03 obtained for a Hamiltonian model. 



VI. SUMMARY 

In this work we found the S2NC-network model to be 
a new model to describe mesoscopic disordered electron 
systems with symplectic symmetry. We constructed the 
topology of the model and the scattering matrices repre- 
senting the potential and spin scatterers. Three indepen- 
dent parameters were needed to characterize the scatter- 
ers. The reflection coefficient r and the transmission co- 
efficient t represent the strength of the spatial disorder by 
defining the mean free path of the network model (in the 
absence of SOI). The coefficient s represents the strength 
of the spin-orbit scattering and defines a corresponding 
spin scattering length. We have shown that our model 
exhibits a localization-delocalization transition. In or- 
der to investigate this transition we calculated renormal- 
ized localization lengths A(r, t, s) by the transfer matrix 
method and obtained the scaling function by an iterative 
fit procedure. The quality of the fit was checked by a x 2 - 
test, which confirmed the assumption of one-parameter 
scaling. We constructed a phase diagram for the system 
showing a metallic phase for any s > 0, if r is small and 
t is large enough. The critical exponent of the correla- 
tion length was obtained to be v = 2.51 ± 0.18. This 
value agrees with previously published results within the 
errors. We pointed out, that in determining the errors 
it is essential to take the correlations into account. The 
critical renormalized localization length was found to be 
A* = 1.83±0.03. By a conformal mapping argument this 
corresponds to a value ao = 2.174 ±0.003 for the scaling 
exponent of the typical local density of states. 
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